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Arising as a fluctuation phenomenon, the equilibrium distribution of meandering steps with mean 
separation (I) on a "tilted" surface can be fruitfully analyzed using results from RMT. The set 
of step configurations in 2D can be mapped onto the world lines of spinless fermions in 1+1D 
using the Calogero-Sutherland model. The strength of the ( "instantaneous" , inverse-square) elastic 
repulsion between steps, in dimensionless form, is [3(/3 — 2)/4. The distribution of spacings s(£) 
between neighboring steps (analogous to the normalized spacings of energy levels) is well described 
by a "generalized" Wigner surmise: pp(0,s) w as 13 exp(— bs 2 ). The value of j3 is taken to best fit 
the data; typically 2 < (3 < 10. The procedure is superior to conventional Gaussian and mean- 
field approaches, and progress is being made on formal justification. Furthermore, the theoretically 
simpler step-step distribution function can be measured and analyzed based on exact results. Formal 
results and applications to experiments on metals and semiconductors are summarized, along with 
open questions, (conference abstract) 
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INTRODUCTORY SYNOPSIS 

For over a decade we have studied the correlation func- 
tions of steps on vicinal surfaces of crystals, that is to 
say crystals that are intentionally misoriented from high- 
symmetry directions. Of special experimental interest is 
the distribution of separations of adjacent steps, corre- 
sponding to ppiO, s) of random matrix theory yj. Tradi- 
tionally, surface scientists have described it by a Gaus- 
sian, with a controversy arising about the relationship 
between the width of the Gaussian and the strength of 
the repulsion between steps. This controversy is resolved, 
and a better representation of the step distribution ob- 
tained, by characterizing the problem with the Calogero- 
Sutherland models and using the simple expression of the 
Wigner surmise, with [seemingly for the first time] gen- 
eral values of (3 in the range of 2 to 10. This range is 
rarely considered by random-matrix theorists. The step- 
step correlation function resulting from the Calogero- 
Sutherland models, for which exact but numerically in- 
tractable solutions exist, can also be applied to data. 

BACKGROUND 

On vicinal crystals there are a sequence of terraces ori- 
ented in the high-symmetry direction and separated by 
steps typically one atomic layer high. If the spacing be- 
tween adjacent steps (the "terrace width") is denoted I, 
then the mean spacing {£) is 1 /tan 0, where </> is the angle 
by which the surface is tilted from the high-symmetry di- 
rection. For theoretical modeling it is convenient to think 
of the crystal as being a simple cubic lattice. (The ex- 
tension to realistic symmetries is not difficult but mud- 



dies the discussion.) By ["Maryland"] convention, the 
steps run on average in the y direction, so that the step 
spacings £ are in the x direction. (The high-symmetry 
terrace plane is then z.) For metals it is a decent ap- 
proximation to take the total energy of the system to be 
proportional to the number of nearest neighbor bonds. 
(For semiconductors corner energies can play a role, but 
these are insignificant for the properties discussed in this 
paper.) Then the ground-state configuration is a set of 
perfectly straight steps (uniformly spaced because of re- 
pulsions we shall discuss shortly). At low temperatures, 
the predominant excitation are kinks in the steps, each 
costing energy e. At higher temperatures, atoms and va- 
cancies appear on the terraces (each with energy 4e in the 
simple model), but these are neglected here. The result- 
ing model is called the terrace-step-kink model (TSK) or 
terrace-ledge-kink model in the literature. 

In this low-temperature picture, the number of steps is 
fixed (by the miscut angle in an experiment and by screw 
periodic boundary conditions in a numerical simulation). 
Thus, the steps never start or stop, and they all have the 
same orientation (say up): there are no anti-steps (or 
down-steps). It is convenient and fruitful to imagine the 
y direction as time-like and to view the configurations of 
steps as the worldlines of a collection of a single kind of 
particle in one spatial dimension (viz. x) evolving in time. 
For physical reasons, the steps cannot cross each other 
(since that would involve atoms suspended above the ter- 
race held with just a couple lateral bonds). Assuming 
also that the steps cannot merge to form double-height 
steps, we can view the evolving ID particles as spinless 
fermions. (Equivalently, in ID, they can be treated as 
hard bosons, a viewpoint exploited in a study of step 
bunching Q. Moreover, the system is highly reminiscent 
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FIG. 1: Illustration of a vicinal generated by a Monte Carlo 
simulation of the terrace-step-kink (TSK) model, showing 
that the elementary thermal excitation is the kink, a unit 
excursion of a step in the x direction. (Adapted from Fig. 
1 of Ref. 27) For experimental STM (scanning tunneling mi- 
croscopy) images of a vicinal surface; see, e.g., Fig. 1 of Ref. 
12, Fig. 1 of Ref. 15, and Fig. 1 of Ref. 34. In the Gruber- 
Mullins approximation, only one step meanders, the others 
being straight along the y direction and uniformly spaced. 
(Cf. Fig. 9 of Ref. 12.) In the fermion picture, the y direction 
is viewed as time-like. In the step continuum model, the steps 
meander continuously rather than by discrete deviations (i.e., 
the kinks are smoothed). 



of soliton lines in 2D incommensurate crystals [3].) 

The obvious next question is how these spinless 
fermions interact. A typical way to characterize the dis- 
tribution of steps is by their terrace-width distribution 
(TWD), i.e. the probability of finding the next step a 
distance £ away. This probability distribution is denoted 
P(s), where s = t/(t) denotes the step spacing normal- 
ized by the [only] characteristic length in the x direc- 
tion. By construction, P(s) is normalized and has unit 
mean. (It ultimately corresponds to p(0, s) of random- 
matrix theory.) For the ground state, P(s) is essentially 
a delta function at s=l. If the steps are imagined as 
straight ( "uncooked spaghetti" ) and deposited randomly 
with probability 1/(1), then P(s) becomes exp(-s). In 
fact, however, at finite temperature the steps do mean- 
der, leading to well-known entropic (or steric) repulsion 
due to the underweighting in the partition function of 
steps on the verge of crossing. This effect is particularly 
noticeable as s approaches 0: P(s) now vanishes with 
power-law behavior rather than growing to unity. For 
s » 1 P(s) must decay faster than exp(-s) to preserve 
the unit mean; analogies with random walkers suggest 



exp(-6s 2 ) as the form of the tail p|. There are many ar- 
guments to show that the entropic repulsion per length 
is proportional to (ksT) 2 / (3£ 2 ] here the step stiffness 
$ = (3(9) + d 2 (3/d9 2 , with [orientation-dependent] free 
energy per length conventionally called (3(9) in surface 
studies. 

In addition to the entropic repulsion, there generally 
is an elastic interaction, also having the inverse-square 
form A/£ 2 : due to repulsion between force dipoles intrin- 
sic to the steps 0-EJ: Viewed from a continuum picture, 
the step involves a pair of sharp angles which could lower 
their energy by healing to form more obtuse angles. This 
process leads to a strain field around the step in which 
atoms in the terraces bounding the step (both above and 
below) tend to move away from the step. Atoms between 
two steps are then frustrated since they are pushed in op- 
posite directions by the pair of steps. Thus, they cannot 
relax fully, leading effectively to a repulsion, with A pro- 
portional to the product of the step dipoles. 

In the fermion picture, the strength of the elastic re- 
pulsion only enters in the dimensionless combination 
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The elastic and entropic repulsions do not simply add 
together: As the strength A of the elastic interaction 
increases, the chance of steps coming near each other de- 
creases, thereby reducing the entropic repulsion. Hence, 
the dimensionless streng^thof the inverse square repulsion 
in essence 71 becomes 
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1 
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In situations where the surface has a metallic surface 
state [one which crosses the Fermi energy] , there can also 
be an oscillatory [in sign] interaction between steps, al- 
beit with a l/£ 2 envelope [HE3- 

Such an interaction 
destroys many of the scaling properties, particularly that 
the form of P(s) is well defined, i.e. that plots of (£)P(£) 
vs. £/(£) are independent of (£). Further discussion of 
this digression is beyond the scope of this paper. 

Knowledge of the elastic strength A is crucial to ad- 
equate characterization of a stepped surface. It is one 
of just three parameters (the stiffness (3 being another) 
of the widely-applied step continuum model [12J and un- 
derlies the "2D pressure" , which determines surface mor- 
phology (e.g., whether steps "bunch") and drives kinetic 
evolution. It is generally believed that the repulsions can 
be adequately approximated as acting between pairs of 
steps at the same value of y, i.e. instantaneously in the 
ID fermion representation. This approximation is clearly 
crucial for viability of the fermion approach. 

The simplest way to estimate P(s) was developed over 
a third of a century ago. In treating polymers in 2D, 
deGennes pointed that the energy of a step (or a polymer) 
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could be described as the line integral over path-length 
ingrements times @(0). Expanding to lowest order, we 
get a constant (the straight step) plus what amounts to 
a kinetic energy proportional to /3 (dx/dy) 2 , the mass-like 
stiffness times the ID velocity squared. 

In the mean-field (i.e. single active s tep with all oth- 
ers fixed at spacings (£)) approximation |l3|, Gruber and 
Mullins (GM) reduced the "free fcrmion" case (.4=0) to 
the familiar elementary quantum mechanics problem of 
solving the Schrodinger equation for a spinless fermion 
of mass oc $ in a hard- wall box of length 2 (£) . The prob- 
ability density of the ground state (the squared ground- 
state wavefunction) is just Pqm-o( s ) — sin 2 (7rs/2), for 
< s < 2; remarkably, the associated ground-state en- 
ergy is precisely the entropic repulsion. For A above 
about 3/2, the problem becomes that of a particle in a 
parabolic potential, i.e. a simple harmonic oscillator, giv- 
ing Pqm oc exp[-(s-l) 2 /2w% M ], with w GM oc i" 1 / 4 ^) 
|14j . These parametric results were dramatically veri- 
fied for Si(lll): not only was the TWD measured at 
one misorientation well fit by a Gaussian, but the TWD 
at a second misorientation was well described by another 
Gaussian whose width was simply rescaled by the change 
in {£}, with no refitting [id llfij. 

However, the Grenoble group pointed out a flaw in 
Gruber-Mullins approach in the limit of very strong A 
[l7j . There the entropic repulsion becomes negligible 
and the steps can be imagined as meandering indepen- 
dently (albeit by a small amount). Then the variance 
of the TWD should not be just Wq M , but twice that 
amount, since both steps bounding a terrace are mean- 
dering. (The actual increase is less than two-fold due 
to lingering anticorrelations.) Hence, a fit of the TWD 
width in the Gruber-Mullins approximation would un- 
derestimate A by a factor of over 3 in the limit of very 
larg e A. Based on roughening theory, the Saclay group 
|18| concluded that the underestimate, for more general 
A, was somewhat smaller (about 2). 

Meanwhile, Ibach noted that the TWD of the free- 
fermion distribution could be far better approximated by 
an expression oc s 2 exp(— bs 2 ), in essence the Wigner sur- 
mise for free fermions, than by Pgm-o(s) = sin 2 (-7rs/2). 
How this comes about becomes evident by recasting the 
preceding model in terms of the celebrated 0] Calogero- 
Sutherland models 



states connected by a non-vanishing matrix elements re- 
pel each other by an amount determined by the symmetry 
of the Hamiltonian. While inapplicable to average quan- 
tities, the approach is appropriate for the fluctuation of a 
large number of energy levels. These fluctuations should 
become independent of the specifics of the level spectrum 
and the weight factors, and so should exhibit a universal 
form depending only on symmetry. This idea can also be 
derived from maximum-entropy arguments [23j |. 

One can draw a correspondence between the fluctua- 
tions of energy spacings between adjacent energy levels in 
nuclei and the fluctuations of spatial separations between 
adjacent fermions in ID systems. (In other chaotic sys- 
tems, as well, there is a correspondence between energy 
and spatial spacings.) The Calogero-Sutherland model 
describes spinless fermions in ID which interact via an 
inverse square potential. Specifically, Sutherland's ver- 
sion |2ij of the Calogero Hamiltonian (2^ for the fluc- 
tuations of uniformly-spaced fermions on an infinite line 
is: 

3=1 3 v 7 l<i<j<JV 

N 

+ - 2 £* 2 , 0) 

in the limits N — > 00 and to — > 0. (In Calogero's original 
model |20|, the last term (oc lu 2 ) is summed over particle 
separations (xj —xA 2 rather than deviations from integer 
positions x 2 .) The ground-state wavefunction for this 
Hamiltonian is 

*o= II K-*/ /2 ex P (^f> 2 ) (4) 

l<i<j<N \ k=l ) 

The ground-state density is recognized as a joint prob- 
ability distribution function from the theory of random 
matrices for Dyson's Gaussian ensembles |24| . 

The Sutherland Hamiltonianj^] similarly describes 
spinless fermions on a circle of radius L (with L — > 00), 
having inverse-square interactions along chords: 



CONNECTION TO RANDOM MATRIX THEORY 
VIA CALOGERO-SUTHERLAND MODELS 

In his study of the spacings of energy levels in nuclei, 
Wigner's starting inspiration was to consider ensembles 
of dynamical systems governed by different Hamiltoni- 
ans with common symmetry properties [i^]. The three 
generic Hamiltonian symmetries are orthogonal, unitary, 
and symplectic. The key ingredient is level repulsion: two 
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n(xj - Xi) 



(5) 

In this case the ground-state wavefunction has the Jas- 
trow form 



* =n 



i<3 



. n(xj - Xi) 



0/2 



(6) 
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With the definition 9{ = 2nxi/L, the ground-state den- 
sity ^>q can be written as 



P(s) 



*8 = IH 



(7) 



i<3 



This ground-state density is again a joint probability dis- 
tribution function from the theory of random matrices, 
now for Dyson's circular ensembles |25| . 

From Eqs. © and JSJ, the dimensionless interaction 
A can be identified as ((3/ 2) ((3/ 2 — 1). Conversely, (3 is 
just the quantity inside the brackets in expression J5|. 
These models can be solved exactly only for the special 
cases (3 = 1, 2, or 4, corresponding to orthogonal, uni- 
tary, or symplectic symmetry of the ensemble. In the 
literature of random matrix theory Q, the TWD P(s) is 
often denoted p(0, s), the probability of finding two lev- 
els separated by s with zero intervening levels. Note also 
that (3 was used in the previous section to denote the 
step free energy per length. (Presumably the letter was 
chosen to emphasis the similarity to but difference from 
the surface free energy per area 7.) Dyson 0,H^| chose (3 
to denote this exponent to indicate the analogy to an in- 
verse temperature (of a Coulomb gas on a circle) . While 
we follow this convention here to minimize the chance 
of confusion, in our applications to stepped surfaces the 
exponent is called g to minimize confusion for that au- 
dience. (The stiffness (3 is essentially unrelated to the 
exponent (3\) 

Random matrix theory leads to exact solutions for the 
ground state properties of the three special cases, in par- 
ticular for P(s). There are even prescriptions that allow 
> gene rate 

dll El- 

in a series and are not useful for fitting experimental data. 
Wigner surmised that P(s) has the simple form 



one to gene rate numerical representations of arbitrary ac- 
curacy THE}. However, these expressions involve terms 



(8) 



for the special cases (3= 1, 2, and 4. The constants ap 
associated with normalization of P{s) and bp producing 
unit mean are: 



a/3 



0+1 



2 



0+2 



and b 







r(fg) 



(9) 



From Eq. JHJ one can readily find analytic expressions 
for the moments and other measurable properties of P(s) 

m 

The argument for Eq. © focusses on the Jacobean 
associated with a change of variables for the Gaussian- 
ensemble probability distribution function p(TL) ~ 
exp[— bN tr(H 2 )] from the eigenenergies to the combina- 
tions of the matrix elements. For an orthogonal ensemble 
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FIG. 2: P(s) vs. s = 
A=0 result (solid curve) 



£/{£} for the exact "free-fermion" , 
the Gruber-Mullins approximation 



sin (ns/2) (long-short dashed curve), and the (3=2 Wigner 
surmise (Eq. (8)) result (dotted curve), barely distinguish- 
able from the exact result). Offset upward by 0.4 for clarity, 
an analogous plot of the exact result for A=2, the Gruber- 
Mullins Gaussian approximation (24/-7T 2 ) 1 / 4 exp(— ^/24(s — 
l) 2 ) (long-short dashed curve), and the (3=4 Wigner surmise 
result (dotted curve). (Adapted from Fig. 1 of Ref. 31). 



((3=1) with N=2, the integrand has a Dirac delta func- 
tion with argument s — [(fan — /122) 2 + 4/i^ 2 ] 1 ^ 2 , which 
vanishes for s—0 only when the two [squared] indepen- 
dent variables do. Hence, P(s) oc s, corresponding to a 
circular shell in parameter space and leading to 0=1. For 
unitary ensembles there is an additional independent pa- 
rameter since h\ 2 becomes (5fte /112) 2 + ($Jm hyi) 2 . Hence, 
P(s) oc s 2 , corresponding to a spherical shell in param- 
eter space, i.e. (3=2. Exact for iV=2, these arguments 
are still fine guides for large N |2^| : As seen most clearly 
by explicit plots as in Haake's text [3(|, P%(s), P2(s), 
and P>i(s) are excellent approximations of the exact re- 
sults for orthogonal, unitary, and symplectic ensembles, 
respectively, and these simple expressions are routinely 
used when confronting experimental data [UIH. (The 
agreement is particularly impressive for i"2(s) and Pa{s), 
which are germane to stepped surfaces.) Furthermore, 
Eq. ||8J has a similar variance at very large (3 to that 
predicted by the Grenoble group, and more generally ap- 
proaches the form of a Gaussian for (3 not too small |3l|. 

The essential innovation in our work is the conjec- 
ture that the generalized Wigner surmise expressed in 
Eq. (JSJ for arbitrary values of (3 provides a good ap- 
proximation for these non-special values. For these gen- 
eral values, there is no recourse to formal justification 
from symmetry arguments. (Indeed, it remains myste- 
rious whether there is some deep underlying reflection 
of the random-matrix symmetry and the corresponding 
Calogero-Sutherland models at the special values of (3.) 
The only way to test our hypothesis is to generate numer- 
ical data. This task has occupied us for some time. We 
have used both Monte Carlo and transfer matrix com- 
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putations of TSK models |3JJ, |32j • The details of these 
studies are beyond the scope of this paper, but the main 
result is that the generalized Wigner surmise does better 
than any of the earlier Gaussian approximations in ac- 
counting for the dependence of the TWD's variance — the 
quantity typically measured by experimentalists — on the 
value of A. In particular, the numerical data confirms 
that the proportionality coefficient between the width w 
of the TWD and A^ 1 / 4 (or more accurately, between w 
and /3 -1 / 2 ) is not constant but increases slowly with A. 

Obtaining the experimentally measurable variance a 2 
of P(s) — essentially w 2 when the Gaussian approxima- 
tion is viable — analytically from Eq. (JHJ), we|3jj can use 
a series expansion to derive an excellent estimate of A 
from a 2 : 



A 



1 

16 



7(a 2 



27 
T 



35 
— c 
6 



(10) 



with all four terms needed to provide a good approxima- 
tion over the full physical range of A, from near zero up 
to around 16, corresponding to /3 ranging from 2 up to 
about 8. The Gaussian methods described earlier essen- 
tially use just the first term of Eq. Ijl0|l and adjust the 
prefactor. 

We have successfully applied these ideas (both direct 
fits to Eq. and use of Eq. (I10|l) to experimental data 
by several different groups [31|, focussing primarily on 
extensive results for vicinal Cu (001) and (111) obtained 
by Giesen, Ibach, and collaborators at FZ-Julich [ssl Isij l . 



PAIR CORRELATION FUNCTIONS 

While the TWD is usually the easiest quantity to mea- 
sure, it is generally not the easiest to calculate since 
it is essentially a many-body correlation function due 
to the demand that there be no fermion between the 
two fermions separated by i (or s). The pair corre- 
lation function his) (also called 1 — Y(s) in random- 
matrix literature[l|) indicates the probability of finding a 
fermion a normalized distance s from a reference fermion, 
regardless of whether there are any between them. For a 
"perfect staircase" this would just be a set of delta func- 
tions at positive integers. More generally, h(s) increases 
like s@ for s <C 1 just like P(s). There are then a se- 
quence of peaks at integer values of s with dips between 
them. The amplitude of the crests and troughs relative 
to the mean density of fermions decreases as s increases. 
For most purposes, but not here, it suffices to use the 
harmonic approximation |35j, I36L l37| (which is essentially 
similar to the Grenoble approximation): 



(i)hp{ S ) = ]T(B m 7T) 1/2 



exp[-B m (s - m) 2 ]; 



(11) 



TT 2 f3 

4B m 



7 + log(27rm) — ci(27rr7i) 



yS-i) j ' +1 (2 ^ w) 

i=i 



2.7 



2j(2j)! 



where 7 ~ 0.577 is Euler's constant and ci is the co- 
sine integral. Note that B m is proportional to /3, so 
that peaks become sharper and higher with increasing 
repulsion. Also, B m decreases slowly but monotonically 
with increasing m. While useful for many applications 
|35l l36l Is?), this approximation for hp{s) proved inade- 
quate for deducing /? from data. 

Forrester [37| found an exact solution for even values 
of p. This expression, which involves Selberg correlation 
integrals, takes nearly a page of text, so is not reproduced 
here. We do note that the ratio of the exact prefactor of 
s@ in hfj{s) to ap in Pp(s) is close to unity, lending further 
support to the utility of the generalized Wigner surmise 
expression 

Subsequently, Ha and Haldane 0, 0, 0] derived an 
exact solution for any rational value of [3. In principle, 
this result also contains information about the "dynamic" 
correlations between fermions (i.e. between steps at dif- 
ferent values of y). Unfortunately, their expression is 
even longer than that of Forrester, and so also is not 
reproduced here. Moreover, it seems to be numerically 
intractable, so of little "practical" use. 

Fortunately, Gangardt and Kamenev recently de- 
rived the following asymptotic expansion for h(s): 



(e)h (s) 

dm = 



00 d 2 (8) 

2 Et^3p« cos (^ s ); (i 2 ) 



it 2 [3s 2 (2tts)« 2 //3 

27rm 

- I II II I I sin I 

(3 



7T/3 



P 



2m 

T 



This formula turns out, remarkably, to be useful 
well below the asymptotic limit, and provides a 
good approximation — much better than the harmonic 
approximation — for s > 1/2; it can be patched onto aps 13 
for s < 1/2. 

In confronting data with Eq. I|13fl . one must beware 
the limited number of steps in experimental images. In 
our application > there were only about a half dozen 
steps per image. Hence, we included an ad hoc linear de- 
cay envelope to account for this limit on h(s) at large s. 
The system in question involved steps on Si(lll) at very 
high temperatures, necessitating a clever device to com- 
pensate for evaporation. For a variety of reasons, the dif- 
ficulty of the measurement led to data that presented ex- 
ceptional challenges for analysis. One problem was that 
steps sometimes disappeared from the image, leading to 
the concern that the TWD might appear wider than it 
in fact was (since sometimes second-neighbor spacings 
rather than adjacent ones might be counted). The inves- 
tigation showed that electromigration from the current 
which heated the sample tended to keep the steps apart, 
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leading to correlation functions indicative of a larger [ef- 
fective] A than predicted from extrapolation of results at 
lower temperature. 

CONCLUDING COMMENTS 

We have seen that steps on misoriented surfaces can 
be represented by fluctuating spinless fermions in 1+1D. 
Thus, step-step correlation functions can be expressed 
in terms of the Calogero-Sutherland model and thereby 
related to random matrix theory. In particular, the ter- 
race width distribution, identified as a joint probability 
distribution, can be well described by the generalized 
Wigner surmise, evidently not only for the special cases 
related to orthogonal, unitary and symplectic ensembles 
but for arbitrary fermion repulsion strength. While per- 
haps reminiscent of some interpolation schemes, such a 
generalization has not, to my knowledge, been made be- 
fore. Furthermore, most of these interpolation schemes 
have involved systems with mixed symmetries (typically 
orthogonal and unitary), with values of (3 between 1 and 
2 (and also between and 1, where /3=0 indicates Pois- 
sonian [random "uncooked spaghetti"] behavior). For 
stepped surfaces the physically interesting range starts 
near (3=2 and goes past 4 up to 8 or even 10. This range 
has received minimal (if any) theoretical attention. 

While the discovery of exact solutions is always intel- 
lectually exciting and captivating, they are often of little 
use in confronting experiments — physical or numerical — 
unless the formulas are numerically tractable. Thus, even 
when a problem is formally solved, neither theoreticians 
or funding agencies should view it as completed until the 
formalism can be rendered in a way that allows the ex- 
traction of numbers, even if approximately. 

A major open question is the "deeper meaning" of 
the generalized Wigner surmise of Eq. JHJ, since there is 
no evident symmetry to argue for the simple expression. 
Howard Richards has taken the lead in trying to derive 
this result from an effective Hamiltonian 44]. This ap- 
proach should lead to the development of the form of 
P(s) for systems with higher-order terms in the repul- 
sion potential and even potentials with an oscillatory 
term such as when metallic surface states mediate the 
interaction. (In systems in which "non-instantaneous" 
interactions (in directions other than x are significant, 
the usefulness of the mapping to ID fermions is likely to 
break down.) It will be interesting to see whether ex- 
pressions similar to the generalized Wigner surmise arise 
in other areas of physics. For example, while consider- 
ing the probability distribution of stock market returns 
in the Heston model with stochastic volatility (in econo- 
physics research), Yakovenko showed that the probabil- 
ity distribution of the variance — which is derived from a 
Fokker-Planck equation, has the form of the generalized 
Wigner surmise 45]! (In this case the value of (3 was 1.6, 



so in a different regime from stepped surfaces. 

While Eq. JSJl assumes that s is continuous, experi- 
ments and numerical simulations involve discrete lattices. 
We have shown 

mi m m 

that so long as (£) is at 
least 4 lattice spacings, there are no significant complica- 
tions due to discreteness and that assuming continuous 
s does not distort the analysis. In particular, the rough- 
ening transition to a facet that occur for discrete lattices 
[Lsl I47I d oes not happen in the range of physical param- 
eters |22j. 

Another skirted issue concerns how many steps in- 
teract. While the entropic interaction ipso facto in- 
volves just adjacent steps, the Calogero-Sutherland mod- 
els jloi HH E3 an d the Saclay approach assume 
energetic interactions between all steps. The Gruber- 
Mullins 0, ^| and Grenoble 0, [23] approximations 
allow all steps or just adjacent ones to interact. Most 
Monte Carlo and transfer matrix simulations assume that 
just adjacent steps interact. In the Gruber-Mullins ap- 
proximation, the curvature of the potential involves the 
inverse 4th power of the distance from neighboring steps, 
so that the effective strength of the repulsion increases 
by C(4) = 7r 4 /90 ~ 1.08 when interactions between all 
steps rather than just adjacent ones are included [3l| : in 
the Grenoble approximation it increases by about 1.10. 
These changes are relatively small, but the good agree- 
ment is curious between the variance predicted by Eq. 
(JHJ and our Monte Carlo results with just adjacent steps 
interacting |32|. 

For surface scientists a fundamental objective is to 
show consistency between the value of A deduced from 
the TWD and that predicted from surface stress using 
Marchenko's formula [||. However, most experimental 
systems do no possess the elastic isotropy assumed by 
that classic result. Hence, this goal has been elusive. 

Another interesting problem is what happens when the 
mean direction of the steps (i.e. y) does not correspond 
to a high-symmetry direction of the surface. We have re- 
cently found that the stiffness computed from standard 
near-neighbor bond lattice models (Ising or SOS mod- 
els) underestimates the experimentally observed stiffness 
by a factor of about 4 for the square-net face of copper 
[46J. The most likely explanations are significant long- 
range interactions or perhaps local relaxations leading to 
3-atom effects. Furthermore, it is by no means obvious 
that the assumption of only "instantaneous" interactions 
between steps is justifiable for substantial azimuthal mis- 
orientation. 
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